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Abstract 

Phase diagrams of the micromaser system are mapped out in terms of the phys- 
ical parameters at hand like the atom cavity transit time r, the atom-photon 
frequency detuning Acj, the number of thermal photons rib and the probability 
a for a pump atom to be in its excited state. Critical fluctuations are stud- 
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values of the detuning we flnd a "twinkling" mode of the micromaser system. 
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1 INTRODUCTION 

There are few systems in physics which exhibit a rich structure of phase transitions 
that can be investigated under clean experimental conditions and which at the same time 
can be studied by exact theoretical methods. Recent developments in cavity quantum 
electrodynamics, and in particular the use of beams of highly excited Rydberg atoms and 
resonant microcavities, has opened up an avenue to study one class of such systems in 
great detail. The micromaser system, which is a remarkable experimental realization of 
the idealised system of a two-level atom interacting with a second-quantised single-mode 
electromagnetic field, provides us with such an example (for reviews and references see 
e.g. Ijl]]). The micro laser [§] is the counterpart in the optical regime. 

Many features of the micromaser are of general interest. It can, for example, be argued 
that the micromaser system is a very simple illustration of the conjectured topological 
origin of second-order phase transitions [^, as will be briefly discussed later in the present 
paper. Various aspects of stochastic resonance can, furthermore, be studied in this system 
The micromaser also illustrates a remarkable feature of non-linear dynamical systems: 
turning on randomness may lead to an increased signal-to-noise ratio |^. 

The basic theory of the micromaser system as developed in Refs. ^ , which is limited 
to the case a = 1 and Au = only, suggests the existence of various phase transitions in 
the large limit. Here is the number of atoms passing the cavity in a single photon 
decay time. A natural order parameter is then the average photon "density" {n)/N, 
where (n) denotes the average photon number with respect to the stationary micromaser 
photon distribution. By making use of an approximative approach in terms of a Fokker- 
Planck equation, details of the various transitions were worked out, like tunnelling times 
between the various phases (see e.g. 0, |[). An exact treatment, in the large A^ limit, 
of the micromaser phase structure and the corresponding critical fluctuations in terms 
of a conventional correlation function of either the atoms leaving the micromaser or the 
photons in the cavity has been given in Refs. , where deviations from the Fokker-Planck 
equation approach were reported. It is of great interest to notice that spontaneous jumps 
in {n)/N and large correlation lengths close to micromaser phase transitions have been 
observed experimentally fl], Most of the experimental and theoretical studies have. 



however, been limited to the case a = 1 and Auj = 0. It is the purpose of the present 
paper to study the phase structure of the micromaser system and extend the results of 
Refs. [y] to the general case with a 7^ 1 and Auj ^ using methods which are exact in the 
large A^ limit. As will be argued below, several new intriguing physical properties of the 
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micromaser system are then unfolded. A short presentation of the main results of this 



paper has appeared elsewhere [|rj . 

The paper is organised as follows. In Section ^ the theoretical framework is outlined, 
and a rapidly convergent large N limit is derived for the stationary micromaser photon 
probability distribution. The phase structure is investigated in Section ^ and the order 
parameter (x) is discussed in more detail in Section ^. The correlation length is analyzed 
numerically in Section |^ as well as in terms of various approximation schemes. The results 
of extensive numerical investigations of trapping states are presented in Section ^. Final 
remarks are given in 0. 



2 THE PROBABILITY DISTRIBUTION 



The pump atoms, which enter the cavity, are assumed to be prepared in an incoherent 
statistical mixture, i.e. the initial density matrix of the atoms is of the diagonal form 

where a + b = 1 . This form of the atomic density matrix is not of the most general form, 
but it leads to the possibility of an exact analytical treatment provided the photon field 
density matrix also is diagonal, which from now on we assume. The rate R of the injected 
atoms is assumed to be high enough to pump up the cavity, i.e. R> ■y = 1/Tcav, or N > 1 
in terms of the dimensionless flux variable = -R/7. Here Tcav is the decay time of the 
cavity and hence 7 is the typical decay rate for photons in the cavity. Furthermore, let r 
denote the atomic transit-time through the resonant photon cavity. 

For our purposes the continuous time formulation of the micromaser system []T^ is 
suitable as a theoretical framework. Each atom is assumed to have the same probability 
Rdt of arriving in an infinitesimal time interval dt, i.e. they are Poisson distributed. 
Provided the interaction with the radiation field of the cavity takes less time than this 
interval, i.e. if r <^ dt, we may consider the atomic transitions to be instantaneous. The 
vector p formed by the diagonal density matrix elements of the photon field then obeys 
the master equation [|12], || 
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where 



L = Lc-N{M -1) . (3) 
Here Lc describes the damping of the cavity (see e.g. Ref. for an excellent account), 



I.e. 



{Lc)nm = {nb + l)[n6n,m — {n + l)6n+l,m ] + nb[ {n + 1)5 

] , (4) 

and M = M(+) + M(— ) has its origin in the Jaynes-Cummings (JC) model 0, ^ where 



M{ + )nm = bqn+lSn+l,m + a{l - qn+l)Sn 



and 



M{-)nm = aqn6n,m+l + b{l - qn)Sn,; 

We observe the unitarity constraints 



(5) 
(6) 



J2 Mnm = 1 , 

n=0 

oo 

^^(-^c)nm = . 



n=0 



The stationary solution of Eq. (|^) is well known 0, and is given by 

nbm + Naqm 



Pn = Po 



Ha 



(1 + rib) 171 + Nbq„ 



where = q{m/N) and 



q{x) 



X 



sin^ (^eVx + A2 



x + A2 

Here we have defined the natural and dimensionless parameter 



where 



(7) 



(9) 



(10) 



6 = ALj/2g 



fill 
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since it is always this particular combination of Acj, and g which is involved, g is 
the single-photon Rabi frequency at zero detuning of the JC model [|1^]. nt is the initial 
mean value of thermal photons in the cavity. Furthermore, the probability is expressed 
in terms of the scaled dimensionless pump parameter 9 = gT\fN . The overall constant 
Po in Eq- (H) is determined by the normalisation condition 

oo 
n=0 

In passing, we observe that at thermal equilibrium, i.e. when 

the probability distribution in Eq. (H) is in fact, as it should be, a Planck distribution. 
This particular distribution depends only on n;,, which is the corresponding mean value 
photon number, i.e. 

oo 

(^) = X] ^ ■ ^^^^ 

n=0 

If Eq. (|T^) is not fulfilled then describes a stationary photon distribution which there- 
fore can be far from thermal equilibrium. 



2.1 Large N Expansion 

The large N behaviour of the stationary photon probability distribution p„ as given 
by Eq. (Q) is difficult to handle. In order to obtain an appropriate expression which is 
rapidly convergent in the large limit, and which neatly expresses the various micromaser 
phases, we can make use of a Poisson summation technique It is then natural to define 
the scaled photon number variable x, defined by x = n/N. The stationary probability 
distribution Eq. then takes the form 

where 



oo 

V{x) = J2 Mx) . (16) 



The effective potential V{x) is expressed in terms of 



Vk{x) = - [ diy\n[w{u)] cos{2TTNkiy) , (17) 
Jo 



where 

rihX + aq{x) 



w{x) 



[1 + Ub) X + b q{x) 

In this expression for w{x), q{x) is given by Eq. (^). We stress that Eq. (|1^) is exact. 

In the large limit Eq. ([T5|) can be substantially simplified since the VA,.(x)-terms in 
Eq. {\r7\j do not contribute in this limit for A; 7^ . It is therefore the nature of the global 
minima of Vq{x) which determines the probability distribution and the micromaser phase 
structure, apart from the zeros of w{x) which correspond to trapping states |16|-||20|. The 
physics of trapping states will be discussed in Section |^. 

2.2 The Thermal Distribution 



If the only global minimum of Vo{x) occurs at x = 0, we can expand the effective 
potential Vq{x) around the origin. A straightforward expansion of Vo{x) then leads to 



Vo{x) = x\n 



l + nb + be^ 



eff 



rib + aOljj: 

where we have defined an effective pump ^-parameter 



2 sin2(M) 



0(x2) 



(19) 



n2 _ n 
^eff - ^ 



(20) 



Eq. (|T9D leads to the following properly normalised thermal photon probablity distribution 

(21) 



Pn 



l + (l-2a) e^f^ 



nb + a 9^ 



eff 



l + Ub + bOj^j \l + nb + b9ifj 



which is convergent provided 

9lff{2a-l)<l . (22) 

Eq. ( |2TD is exact in the large N limit, and corresponds to an increase of the temperature 
in the cavity. If a < 1/2 + the distribution Eq. (pT]) is therefore always valid. If 
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A = convergence requires that (P' (2a — 1) < 1, i.e. Q has to be sufficiently small when 
a > 1/2. We find that Eq. ( pTj ) also can be used to compute expectation values for 
finite A^, even though the accuracy then depends on the actual micromaser parameters. 
When the micromaser is described by the thermal distribution in Eq. (0), we say that 
the system is in the thermal phase. 



2.3 The Gaussian Distribution 

If, on the other hand, non-trivial saddle-points of Vq{x) exist, Eq. (|15D can be written 
in the form 

Pi^) = ^Pji^) , (23) 

j 

where Y2j ^ ^^^^ saddle-points corresponding to minima of Vq{x), and where 
pj{x) is p{x) for X close to the saddle-points x = Xj. The saddle-points, x{9), are deter- 
mined by Vq{x) = 0, since, as noted above, the Vfc(x)-terms in Eq. (|T^) for A; 7^ do not 
contribute in the large iV-limit. For a fixed a and A we then have 

(2a - 1) sin' (eV^d) + ) - i^i^) + A') = , (24) 

corresponding to w{x{6)) = 1. This equation is independent of rib and has non-trivial 
solutions only when a > 1/2 + A'/2. It is clear that a, A and x{6) at a saddle-point are 
restricted by the condition 

x{9) + A' < 2a - 1 , (25) 

for any 6. The saddle-points x = x{9) can conveniently be parametrically represented in 
the form Q 

x(0) + A' = (2a - 1) sin' , (26) 



g(0)= , ^ , (27) 

V2a-1 |sm0| ^ ^ 



with 



> 00 = arcsin(|A|/V2a - 1) . (28) 



We define branches of 0, labelled by k = 0, 1, 2, ... such that varies in the intervals 



(po + kn < (p < (k + 1)71 - (po , A; = 0,1, 2,... . (29) 

Except for the first branch k = 0, the saddle-points corresponding to each of these 
branches are double- valued, that is, there are at most two values of x{6) for a given 
6. The first branch = is single-valued. 

For a given a, A and 6, let saddle-point corresponding to a 

minimum of Vq{x). The equilibrium distribution in the neighbourhood of this minimum 
value Xj is then given in the following Gaussian form 



Pj[x) 



'^L= g-f ^o"(^.)(^-^.)' . (30) 



Here Tj is determined by the normalisation condition for p{x), i.e 



and where 



p-NVo{xj) 

J". = 1 (31) 



T/'V ^ - (2^-1)^ g(x) - xq'jx) 

^ol^J-^ + ^^(2a-l) x2 • ^^^^ 



For the given parameters, the sum in Eq. ( pT]) is supposed to be taken over all saddle- 
points corresponding to a minimum of Vo(x), i.e. all saddle-points corresponding to 
Vq{x) > 0. If a; = Xj does not correspond to a global minimum for the effective po- 
tential Vo{x), then Tj is exponentially small in the large limit. If a; = xj corresponds 
to one and only one global minimum, we can neglect all the terms in the sum in Tj but 
m = j, in which case Tj is reduced to Tj = ^/\^Jxj). In the neighbourhood of such a 
global minimum Xj the probability distribution in Eq. (|30|) is therefore reduced to 



which is, strictly speaking, only valid in a sufficiently small neighbourhood of Xj. 

However, since Pj{x) in Eq. (|30|) is exponentially small when x is not in the neigh- 
bourhood of Xj, the probability distribution for any x is given by 

Pix) = ^Pji^) > (34) 
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where Ylj* denotes the sum over the global minima of Vq{x) only. Hence, when there 
are several saddle-points, the actual maser phase is described by the saddle-points which 
correspond to the global minima of Vo(x) only. In the next section, we will actually argue 
that, for a fixed a, A and 6, we can have at most two competing global minima in the 
large limit, i.e. there can be at most two terms in the ^^.-sum in Eq. (^4]). 

When the micromaser is described by a distribution like Eq. (^) the mean occupation 
number (n) is proportional to the pumping rate N (see Section ^). The cavity then acts 
as a maser, i.e. the system is in a maser phase. 



3 THE PHASE DIAGRAM 



The probability distribution Eq. (^) determines micromaser phases as a function of 
the physical parameters at hand. In this section we will map out phase diagrams in the 
a- and 6'-parameter space for a given 7^ and A. In general a phase diagram is then 
determined by mapping out the global minima of the effective potential Vo{x). 

By the substitution = + A^, the effective potential in Eq. (0) can be written 
in the form 

Vo{x) = Vo{(P,e) = -^ r dxxHwix,0)] . (35) 
^ Je\A\ 

Here we have defined 

w{x,0) = — , (36) 

l + nb + bq{x,0) 

and 

q{x,e) = e^'^ . (37) 
X 

In Eq. (^), the upper integration limit and the pump parameter 6 may take on arbitrary 
values. However, by choosing x(0) and 6'(0) according to Eq. ( p6D and Eq. ( |27| ) the 
effective potential is always at an extremum. We then arrive at an effective multi-branched 
potential as given by Vo(0) = Vo(0, 0{(p))- For a given branch k, the variable (j) is limited 
by Eq. m). 
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If we regard Vo(0) as a function of 9, each branch of Vq{9) = Vo(0, 6'(0)) is at most 
double- valued except for branch k = 0. One sub-branch then corresponds to a maximum 
(V^"(a;(0)) < 0), which is swept out first as increases, and the other corresponds to a 
minimum (V^"(x(0)) > 0). Here the Vo"(^(0)) is explicitly given by 

1 — d) cot (b 

W0))= . 2 .r / J .,1 • (38) 

sm [a + nh[2a — IjJ 

To reconstruct the actual value of Vo(x) for the A;:th minimum (maximum), we write 
Vo(x) = Vo(0), evaluated for the sub-branch with Ko"(x(0)) > (V^o"(a;(0)) < 0). 
A convenient representation of Vo(0) for the branch k (given a, rib, A) is 

Vo(0) = [a + nfc(2a-l)]J(0) , (39) 

where 



r-cj) 



dxx' sin(2x)-2sin (x)/x _ ^^q^ 

^o+fcTT a6'(0)2sin^x] [x^(l + n^) 66^(0)2 sin^ x] 



Here the upper integration limit has to be chosen according to Eq. (^). We notice the 
characteristic pre-factor a + nh{2a — 1) in Eq. (^Uj) which actually determines the typical 
dependence of the parameters a and rif, at a saddle-point. 



3.1 The Critical Parameters and 9k 

As a function of 9, the first extremum of Vo(x) occurs at 9q = 6'(0o), where 0o is given 
by Eq. (|28|). The critical pump parameter = 0Q(a, A) is therefore given by 



n* arcsin(|A|/v/2a- 1) 

° = |X| ' (^^) 



which is equivalent to 



^ = 7^ + 77—2 77^ ■ (42) 



1 A 

2 ^ 2^h?(^*A) 

When the system is not detuned, i.e. when A = 0, this equation reduces to 
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Eq. ( ^Tj ) determines the first thermal-maser critical line in the a- and 9- phase diagram. 
In the maser phase, the order parameter (x) = {n)/N approaches zero when 9 approaches 
9l{a,A). Furthermore, in the large N limit, (x) is always zero in the thermal phase 
(see Section ^). Hence, the order parameter is continuous on this critical line. The first 
derivative of (x) with respect to 9 is, however, discontinuous. In the maser phase we 
actually have 

^|>=2(^2^f-ii^ . (44) 

d9 tan0 — cp 

with > 00- When approaching the critical line from the maser phase, we see from 
Eq. ( ^Ij) that d{x)/d9 is non-zero (except for a = 1/2). The first thermal-maser critical 
line therefore corresponds to a line of second-order phase transitions. 

In passing, we mention that in a topological analysis of the second-order phase tran- 
sitions of the micromaser system, Vo{x) will play the role of a Morse function. The 
configuration space Ai is then the one-dimensional space of x{= n/N). The sub-manifold 
Ai^ defined by 

My = {x e M \ p{x) ~ exp( -NVq{x) ) <v} , (45) 

then describes the change of topology close to the second-order phase transitions 9^. This 
will be true for all second-order transitions discussed below. We observe that for finite 
N the change in topology, i.e. the appearance of an additional disjoint x-interval due 
to the appearance of a new local maximum in p(x), in general occurs before the actual 
second-order phase transition. 

Critical points where new extrema of Vo(x) appear are determined by Vq"(x) = 0, i.e. 
non-trivial solutions of 

tan = , (46) 

independent of the physical parameters of the micromaser. This equation has infinitely 
many positive solutions = 0jt, where k = 1, 2, ... . Corresponding to these solutions we 
have the critical pump parameters 9k = 9{(f)k), i-e 

9k = J—, -^ , A; = 1,2,3,... , (47) 



^/2a — 1 I sin < 



>k 



independent of and A, for which the fc'th branch comes into existence. For a given 
branch k > 1, the V"q"(x(0)) > sub-branches and the VJ)"(x(0)) < sub-branch coin- 
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Figure 1: The order parameter (x) = {n)/N, as a function of 9, when a = 1, 
A = 0, nft = 0.15 and AT = 1000. A dotted vertical line at 6* = 6'fe, where 6'o = 6*0' 
indicates where the k:th branch comes into existence. 9q corresponds to a second- 
order phase transition. The numerical values of the pump parameters are: 6^ = 1, 
9i = 4.603, 92 = 7.790, 9^ = 10.950, ^4 = 14.102 and ^5 = 17.250. The various 
branches at the extremum of the effective potential, i.e. Vq{9) = 1^(0, are 
also shown. The intersection between two neighbouring VJj"(x(0)) > sub-branches 
at ^ = Oli^_^i is marked by a circle when this coincidence occurs at a global minimum 
of Vo{9). At these first-order maser-to-maser phase transitions the order parameter 
(x) makes a discontinuous change in the large N limit. The numerical values of the 
transition parameters , 1 are: ^qi ~ 6.6610, 0*2 ~ 12.035 and ^ 17.413. 
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Figure 2: The order parameter (x) = {n)/N, as a function of 9, when |A| = 0.5 
and ah other parameters as in Fig. |l[ A dotted vertical line at 9 = 9j^, where 
^0 = ^0' indicates where the A;:th branch comes into existence. The numerical value 
of 9q is 9q = 1.047. The numerical value of the 0^, A; 7^ 0, is the same as in Fig. || 
since this critical 9 only depends on a. The various branches at the extremum 
of the effective potential, i.e. Vq{9) = VQ[(j),9{(f)))^ are also shown. First-order 
thermal-to-maser phase transitions occur at the pump parameters 9^^^ and second- 
order maser-to-thermal transitions occur at 9q, 9q^ and ^^j. The circle at ^ = 
in the figure indicates a maser-to-maser transition. The numerical values of the 
critical transition parameters in the figure are: ~ 1.047, 9q^ « 5.236, 9^1 « 6.193 
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£2 



11.906 and 6*, 



23 



17.425. 
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cide at this particular critical point 9k (see e.g. Fig. The critical point 9 = Ok also 
corresponds to x'{9k) = oo as seen from Eq. (0). 

Each sub-branch with Vo'(a;(0)) > of the effective potential Vo{9) has one and only 
one minimum. This can be seen from the expression 



dVo{9) 
d9 



[a + n5(2a-l)] J(0) , (48) 



where 



and Vq{9) = Vo(0, 9{(f))). The upper integration limit must then be chosen according to 
Eq. (^9|) where we, of course, only consider the V^"(x(0)) > sub-branch for k > 1. 

We can now study the actual intersections of the various branches at a common global 
minimum by making use of numerical and analytical methods. The procedure is, for 
example, illustrated in Fig. |I] and Fig. ^. We find that Vo(^) is such that if two branches 
do intersect at a common global minimum, these branches correspond to k and k + 1, i.e. 
they are consecutive branches. A given branch can also intersect with Vo(^) = 0, i.e. the 
thermal phase. In addition two consecutive branches can intersect when Vo(^) = 0. We 
then have a triple point. 

3.2 The Critical Parameters Olj^^^, O^f^ and 



If, for a given a, rih and A, we have a transition from the maser branch k to the 
neighbouring maser branch k + 1, then there exist a (pl and such that 

= K)(0^+i) and 9{(I>1) = 9{<I>1^,) ^ 91,^, , (50) 

where 

^oi<Pl) = -J7^r d4><p\n[w{cP,9{cPl))] , (51) 

and where (pl is in the interval 0^ < 0^ < (tt — 0o) + kn. The reason why we chose 
(pl > (pk is that we are looking for solutions of Eq. (^) for which Iq"(x(0)) > 0. The 
pump parameter according to Eq. (|50D can be expressed as 
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I sin 0*1 • ^^^^ 
Furthermore, for this particular combination of the parameters, the order parameter (x) = 
(n) /N makes a discontinous change as can be seen from the Eq. (p6|), i.e. the discontinuity 
in (x) is given A{x) = x{(f)l_^-^^) — x{(f)l). We therefore have a first-order phase transition. 

If, on the other hand, a given maser branch corresponding to a global minima does not 
intersect with a neighbouring maser branch, it is the intersection with the thermal branch 
which determines the phase transition (see e.g. Fig. For a given branch k > 1, let 
6^1^ = 0{(f)lf}) denote the pump parameter of this thermal-to-maser transition. The value 
of 

4'*tk is then determined by the solution of 

K)(0aO = O , (53) 

where 

M<Ptk) = -7ui^ d<p<PHw{<p,e{<p;,))] . (54) 

^mk) Je{<p*^)\A\ 

We determine numerically, where 0^ < <Ptk < (tt — 0o) + kn and k = 1, 2, 3, ... . The 
corresponding pump parameters are 6*^^ = 6'(0^^), i.e 

^r. = ^^n%M ' ^ = l'2,3,... . (55) 
V2a - 1 |sm0^fc| 

Due to the same reasons as discussed above, such thermal-to-maser transitions correspond 
to first-order phase transitions. 

Furthermore, let $1^ = 9{(f)\^ denote the maser-to-thermal transition for a given branch 
/c > (see e.g. Fig. |). The value of is determined by one of the solutions of 



^o(0h) = O , (56) 



where 



M^lt) = -7uJrY2 #01n[i^(0,^(0*J)] . (57) 

kt) Je{<t>it)\^\ 



We see that 



'kt 



(vr - 0o) + k-n , (5J 



15 




e 



Figure 3: The phase diagram for the micromaser system when nf, = 0.15 and 
A = 0. All the critical lines converge to 1/2 in the large 9 limit. The thermal-to- 
maser critical line OQ{a) is determined analytically by Eq. (^). The other critical 
lines which are maser-to-maser transitions, are determined by Eq. (|K 



for all branches k > 0. When A 7^ 0, Eq. (|56D is trivially fulfilled since the upper 
integration limit is equal to the lower. If, on the other hand, A = 0, then Eq. (^) is 
fulfilled since 0{(f)l^) = 00. The pump parameter 91^ = 0{(f)kt) is given by 



|A| 



This equation is equivalent to 



1 A^ 

a = - + ^- , A: = 0,1, 2,... , (60) 

2 2sin2(0*^A) > > > > K J 

which is of the same form as Eq. (^). Due to the same reasons as discussed above, such 
a maser-to-thermal transition corresponds to a second-order phase transition. 

Equipped with these results, we can now construct a complete phase diagram in e.g. 
the a- and 6'-parameter space for a given and A. 
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3.3 Phase Diagram 



The phase diagram when = 0.15 and A = is shown in Fig. ^ The first critical hne 
in this figure is given by Eq. (^). As aheady mentioned, this corresponds to a second- 
order thermal-to-maser phase transition. In the region above this hne the mean occupation 
number (n) grows proportionally with the pumping rate A^. The cavity therefore behaves 
like a maser in this regime. For values of a and 6 below this line the only global minimum of 
Vq{x) occurs at X = 0. In this particular region the probability p„ is given by the thermal 
distribution in Eq. i.e. the micromaser is in the thermal phase. We remark too that 
Eq. ( ^2]) also determines the radius of converegence of the thermal probability distribution 
in Eq. (0). The other critical lines in Fig. ^ are determined by Eq. (pOD since for rih = 0.15 
and A = 0, all neighbouring VQ"(a;(0)) > sub-branches of Vo{9) = Vo{(f){9),9) intersect 
for some a > 1/2, see Fig. |l]. 

When the micromaser is detuned the phase diagram is more complicated. Fig. § shows 
a typical example of a phase diagram when A 7^ 0. As we can see from this figure, the 
first two maser phases are well separated. The critical value of A for such a separation 
of phases is in general determined by considering phase separation on the line a = 1. 
For a given n^, let A = Afe^+i be the corresponding critical value for phase separation. 
With the definition = 0o, where 0o is given in Eq. (0), A^^+i is determined by the 
transcendental equation 

^(0 = Wfc) , fc = 0,l,2,... , (61) 

i.e. the solution of 

|A| = ^^^^[(A; + l)7r-arcsin|A|] . (62) 

Here 0^^, is determined numerically according to Eq. (|55|). For Ub = 0.15 the first critical 
value of detuning is |Aoi| ~ 0.408. When the detuning is larger than Akk+i, branch k will 
separate from branch k + 1 (see e.g. Fig. |^). In passing, we notice that A^k+i is restricted 
by Alj^_^_-^ < 1, for all possible branches. 

The first critical line in Fig. ^ corresponds to a second-order thermal-to-maser (maser- 
to-thermal) transition to the left (right) of its minimum. The minima of the critical lines 
are marked by short vertical lines. For any other critical line in Fig. ^ the transition is a 
first-order thermal-to-maser (second-order maser-to-thermal) phase transition to the left 
(right) of its minimum unless it intersects with another critical line. In the region between 
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a 0.7 




e 

Figure 4: The phase diagram for the micromaser system when Uf, = 0.15 and 
|A| = 0.5. In this diagram we have plotted the range of vahdity of the thermal 
distribution (dotted lines). When not visible these dotted lines overlap with the 
solid critical lines. The minima of the critical lines are determined by the condition 
sin^(A0) = 1 and are marked by a short vertical lines. The first critical line OQ{a) in 
this phase diagram is determined analytically by Eq. The critical lines 0^i^{a) 

are given by Eq. ( [55| ) and Oltia) are given by Eq. (59). Triple points are indicated by 
circles. The numerical values of these triple points are (a23i ^23)*^*^'*^ = (0.98, 17.78) 
and (034,^34)*"^'^ = (0.96,24.04). 



any thermal-to-maser line and the neighbouring dotted line, the thermal probability dis- 
tribution does not correspond to a global minimum of Vo{x). The phase diagram in Fig. ^ 
also contains triple points. These points are marked by circles. They are mathematically 
determined by 

Voi<j);,) = Vo{<l>l,) = and ^(0,*,) = e(</.*J . (63) 
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4 THE ORDER PARAMETER 



When the system is in a maser phase the probabihty distribution is given by Eq. (pBI). 
(x) = {n)/N is then, of course, trivial to compute. For a given a, A and 6 (x) is simply 
given by 

{x)=x, , (64) 

where xj is the value of x where Vo{x) reaches its global minimum, i.e. xj is one of 
the saddle-points as determined by Eq. (P^. In the maser regime, the photon number 
behaviour (see e.g. Fig. |I|) is well-established in the literature (see e.g. Refs.|I[]). 

For values of a and 6 corresponding to the first critical line Eq. (^), the order param- 
eter (x) is zero. Mathematically this follows trivially by substituting 0o in Eq. (^) into 
Eq. ( p6[ ) and using Eq. (Q). Below this critical line, the micromaser is in the thermal 
phase and the order parameter is then given by 



■ni, + a9 



which is zero in the large limit. We realize from Eq. (^) that the photon number (n) 
is in general a periodic function of 6 when the system is detuned. The period of (n) is 
then |A|6'. It reaches its maximum (n^ + a/A^)/{l + (1 - 2a) /A'^) for A9 = {n + l/2)7r, 
where n = 0, 1, ... . The minimum of (n) is Ub and occurs for |A|^ = nvr. In passing, we 
note that the maximum value {rib + c^/^^)/(l + (1 ~ 2a)/A^) is reduced to in thermal 
equilibrium, for any Ub, A and 6. We also note that when nb = the cavity is dark, i.e. 
(x) = for every integer multiple of 9 = vr/|A|. 

When the system is not detuned, Eq. (|55|) is not oscillating and it then reduces to 



= l + (l-2a)g^ • ^''^ 
In the large 9 limit, Eq. (|66|) converges to a/(l — 2a) (see Fig. H). 

The order parameter (x) can also be studied as a function of the probability a for 
the pump atoms to be in an excited state. Varying a over all allowed values keeping 
other physical parameters fixed, the order parameter (x) makes a discontinuous change 
for every combination of the physical parameters corresponding to a critical line. Hence, 
(x) as a function of a exhibits a plateau-like behaviour when N is sufficiently large. Such 



19 




Figure 5: This figure shows the different behaviours of Eq. (pSl). rih = 0.15 for 
ah curves. The straight line corresponds to a = n{,/(l + 2nh) = 3/26, i.e. thermal 
equilibrium. The sohd curve corresponds to a = 0.2 and |A| = 0.5. The fat dotted 
curve corresponds to a = 0.2 and A = 0. The dashed-dotted curve corresponds to 
= and A = 0. 



a behaviour is illustrated in Fig. ^. In this figure we observe a rather slow convergence 
to this plateau-like behaviour. 

4.1 Twinkling and Detuning 



Above we have observed that the mean value (n) may oscillate as a function of 6 when 
the micromaser is in the thermal phase. We say that the system is then in a twinkling 
mode. This twinkling behaviour has a close resemblance to the observed atomic inversion 



revivals 21, 22 



A somewhat different twinkling behaviour of the micromaser system occurs when the 
thermal phase intersects with a maser phase. This feature is illustrated in Fig. ^ with 
|A| = 0.5 and n?, = 0.15. For a given a > 1/2 + A^/2 we then see that the system has 
repeated thermal-to-maser and maser-to-thermal transitions as 6 increases. This twinkling 
behaviour is not strictly periodic in 6. The winkling phenomena will now, however, be 
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Figure 6: The order parameter (x) = {n)/N as a function of a when rif, = 0.15, 
A = 0, = 25 and N = 100, 500, 1000. Each step in (x) corresponds to a point on 
one of the transition curves in Fig. 0. 



more pronounced since, for large N, the maser will be "dark" in the thermal phase, but 
(n) is large in the maser phase since (n) is proportional to N in this region. 



5 THE CORRELATION LENGTH 

In this section we study long-time correlations in the large limit. This notion of 
a correlation length for the micromaser system was first introduced in Refs. 0. These 
correlations have a surprisingly rich structure and reflect global properties of the station- 
ary photon distribution. The probability V{s) of finding an atom in a state s = ± after 
the interaction with the cavity, where + represents the excited state and — represents the 
ground state, can be expressed in the following matrix form : 

r{s) = u^'^M{s)p , (67) 
such that P(+) + ) = 1. The elements of the vector p are given the equilibrium 
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distribution in Eq. (H), and the matrix M(s) is given by Eqs.(||) and (H). The quantity 
is a vector with all entries equal to 1, = 1. Furthermore, the joint probability for 
observing two atoms in the states si and S2 with a time-delay t between them, is given 
by i 

= M°^M(s2)e-^^*5(si) p , (68) 
where L is given by Eq. (|^) and where 

S{s) = {l + Lc/N)-'M{s) . (69) 

We observe that V{+,—,t) = V{—,+,t) P, which means that the cavity photons and 
the outgoing atoms are not quantum-mechanically entangled. Statistical correlations do, 
however, exist. 

A properly normalised correlation function is then defined by 

7 it) = , (70) 

where {ss)t = J2si S2 ^i^^'Pisi, S2,t) and (s) = J2s sV{s). This correlation function satis- 
fies — 1 < 7'^(t) < 1. At large times t ^ oo, we then define the atomic beam correlation 
length U by i 

7^(t)~e-*/?^ . (71) 

The lowest eigenvalue A = of L determines the stationary equilibrium solution p. The 
next non-zero eigenvalue Xm of L, on the other hand, determines the typical time scales for 
the approach to the stationary situation. This eigenvalue can be determined numerically. 
The relation between A„2 and the atomic beam correlation length is 1/E,a = l^m- For 
photons we define a similar correlation length ^c- It follows that the correlation lengths 
are identical, i.e. = = ^ 0- 

The correlation length 7^ is shown in Fig. ^ (A = 0) and Fig. |] (| A| = 0.5) for various 
values of A^. Furthermore, Fig. ^ shows 7^ for various values of the detuning. When the 
detuning is sufficiently small, we observe from these figures that the correlation length 
exhibits large peaks for different values of 6. In the large A^ limit, numerical studies reveal 
that these large peaks occur at the transition parameters 9q, 9lf^^-^ and/or 6*^, depending 
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Figure 7: The logarithm of the correlation length 7^ as a function of 6 for various 
values of = 25, 50, 125 when n?, = 0.15, A = and a = 1. The numerical values 
of the vertical lines are as in Fig. ||. 

on the values of the given and A (see e.g. Fig. ^ and Fig. ||). We will discuss the 
behaviour of these peaks in more detail in Section ^.4| . When the detuning is sufficiently 
large, i.e. > 2a — 1, the correlation becomes smaller and behaves in a strictly periodic 
manner as a function of the pump parameter 9 (see e.g. Fig. 

In order to arrive at a better understanding of the behaviour of the correlation length, 
we will derive various approximative expressions in the following subsections. We will 
use three different methods in order to derive such expressions for ^ and we will compare 
them with the exact numerical results. 

5.1 The Eigenvalue Value Method 

By making use of the orthonormality conditions for the left (m„) and right (p„) eigen- 
vectors corresponding to the eigenvalue problem as defined by Eq. (H), i.e. uFL = Xu^ 
and Lp = Xp, it is shown in ^ that an eigenvalue satifies the equation 
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A = ^ PnBn{Un - U^-lf , (72) 
n=0 

where 

Bn={n, + l)n + Nhqn • (73) 

and Pn = P^n^n- The left eigenfunction Un corresponding to the first non-zero eigenvalue 
\nz has one node. It is natural to assume that this node occurs at (n). This is so since in 
the neighbourhood of (n) we can then use the Ansatz 

n — in) , , 

Un= — , (74) 

where a\ = {ri^) — (n)^ as usual. The average values are assumed to be taken over the 
equilibrium distribution p„. It then follows that u-p = {un) = and u-p = {uD = 1. By 
the Ansatz Eq.(^) it now follows that 

2 

7^ ^ 7ei,(e) = -— — — -- , (75) 
{Ub + l){n} + Nb {qn) 
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since 7^ = 1/A„2- Eq. (^) can be simplified in the large limit. In this case we have 
~ 9ljj-{n)/N. If, in addition, the micromaser is in the thermal phase, then Eq. ([75| ) 
reduces to 

which explains the periodic behaviour of the correlation as mentioned in the beginning 
of this section. We also observe that 7^^; = 1 when a = 1/2, independent of any of the 
other physical parameters at hand. 



5.2 The Master Equation Method 

Let us now derive an approximative expression for the correlation length by making use 
of another method. For a general right eigenvector of the matrix L we define p{x) = Npn 
and write it as p{x) = p{x)u{x) with the left eigenvector u{x) = Un- It is then shown in 
g that 
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A u{x) = [x — (2a — 1) q{x) ] 



du{x) 
dx 



1 d 

Ndx 



ribX + aq{x) 



^ du{x) 
dx 



(77) 



in the large limit. Here we again make use of the fact that the eigenfunction corre- 
sponding to Xnz has only one node, say aX x = xq. The function u{x) is therefore given 
by u{x) ^ X — Xq in the neighbourhood of Xq. Furthermore, when 9^yx{9) + ^ 1, we 
can expand q{x) around x = 0, in which case Eq. ([77|) reduces to 



2 liO) 



N 



where 



a{e) = 1 - (2a - 1) 9' 



eff ' 



P{9) = (2a- 1) 



sm^{9A) 0sin(0|A|)cos(0|A|) 



A^ 



l{9) = nb + a9lff 



A3 



(78) 

(79) 
(80) 
(81) 



Notice that Eq. (^) is only valid when we are not too close to a critical line. For small 
perturbations around xq, i.e. x = xo + 5x, the eigenvalue Xnz is then 



\nz = a{9) + 2P{9)xo , 
where Xq is easily solved by the aid of this equation and Eq. (|78D : 



(82) 



x,{9) 



a{9) 



a 



id) Y i{9) 

2m) mN 



^3) 



Since 1/^ = 7A„2, we get 



a\9) + 



N 



In the large limit, where Xq = (see Eq. (|83D ), Eq. (|8^) is reduced to 7^a//(6') 
i.e. 



(84) 

l/a{9), 



liM{9) 



1- (2a -1)02^^ ' 



^5) 



in agreement with Eq. 
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Figure 10: Comparison of log(7^) as a function of 6 for various analytical 
expressions of the correlation. The parameters for all lines are a = 1, |A| = 0.5, 
Ub = 0.15 and = 100. The solid curve is the exact correlation length. 



The correlation length iuid) exhibits its first peak at = (see e.g. Fig. [T^ and 



Fig. |TTD. The correlation length Eq. (p^ at this particular peak is 



2a- 1 



where 



N 



a + nf,(2a — 1) 



7) 



1 — X cot X 

In Eq. (^) the probability a can only be chosen in the interval l/2 + A^/2 < a < 1. 
This equation shows that the correlation 9 = 6^ grows as \/N. Again we notice the factor 
a + nf){2a — 1). Also note that Eq. (^) vanishes when a = 1/2. Since by assumption 
e\A\ < 1, thenEq. (|7D isx(^|A|) ~ V^[l-{eAf/5]. The correlation length in Eq. (|8|) 
is then in agreement with the results of Ref.[p| for A = 0. 
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Figure 11: Comparison of log(7i^) as a function of Q for various analytical 
expressions derived in the text. The parameters for all curves are a = 0.75, |A| = 
0.5, nft = 0.15 and = 100. The solid curve is the exact correlation length. 



5.3 The Mean Field Approximation 



We can also use a mean field approximation in order to get an approximative and an- 
alytical expression for the correlation length ^. When the atoms have Poisson distributed 
arrival times, the continuous master equation gives the following exact equation for the 
average photon occupation number 

+ a(g„+i)-%„) . (88) 

By making use of the mean field approximation ^ qi^n) = qH^)), the above equation 
simplifies considerably. The stationary solution, (x) = Xq, of Eq. (|88| ) is determined by 
a transcendent equation which, generally, can be solved numerically only. However, with 
the approximation 

q{{x)) ^ q{{x) + 1/N) ^ q{{x) + f/N) , (89) 



i a\x) 



dt 
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Figure 12: Comparison of log (7^) as a function of 6 for various analytical 
expressions of the correlation as derived in the text. The parameters for all curves 
are a = 0.25, |A| = 0.5, Uh = 0.15 and N = 100. The solid curve is the exact 
correlation length. In the large N limit all curves overlap completely. 



the stationary solution xq has a simple parametric representation. To get a consistent 
expression for the correlation with Eq. (^) the weight factor / has to be 



/ 



With the approximation Eq. 
to 

1 d{x) _ 
7 dt 



2a - 1 

(|89|), which clearly is very good for large A^, Eq. 



{x)-^_+{2a-l)q{{x) + f/N) 



(90) 
reduces 

(91) 



The stationary solution xo of Eq. ( pi]) then satisfies the transcendental mean field equation 



Xo 



^ + (2a-l)g(xo + //iV) 



(92) 



In the large limit, Eq. (p^ ) reduces to Eq. (^) after a multiplication by Xq. In contrast 
to the saddle-point equation Eq. (|^), Eq. (p2D therefore always has the trivial xq = 
solution in the large limit. The general solution of Eq. (l9^ ) can be written in the 
following parametric form: 
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where 



(93) 



h{(f)) = L_pL + A' - (2a - 1) sin^ , (94) 
^?(0)-^ + A^^ + ^(2a-l)sm^0 , (95) 



and 



g(0)= , ^ . (96) 

v/a;o(0) + //iV + A2 

For small perturbations around Xq, i.e. (x) = Xq + we find the equation of motion 

-[l-(2G-l)g'(xo + //iV)] 5x , (97) 



7 c?t 
where 



- ^ - (^fs^ .„=(9vOT) (98) 



+ 7 ^ sin(^Vx + A2) cos(0Vx + A2) . 

(x + A^j'''"^ 

In terms of 0, Eq. (^) at x = Xo(0) reads 

g^('^)= (^^(0) + //iV + A^)^ [A^in^0 + xo(0)0sin0cos0] . (99) 
From Eq. (^7p we immediately see that the correlation is given by 

7e^7eMF(0) = . . (100) 

1 - (2a- l)g'(0) 

Notice that 7^mf is here a function of the parameter 0, in contrast to Eq. (|75|) and 
Eq. (|§. 

When the micromaser is not detuned, Eq. ( |100| ) reduces to 

76/f(0) = — M • (101) 

^ ' 0sm0cos0 

1 — 2a — 1) 7. TF 

^ ^ (2a - 1) sin^ + ^ 
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The peak at 6 = 6q corresponds to 

3K + /)\'/' 



^iV(2a- 1) 

Substituting Eq. (|102|) into equation Eq. ( [101| ) the correlation is 



(102) 



. *^ 2a - 1 3iV , , 

2 y a + nb[2a — 1) 

for some a the interval 1/2 < a < 1. This correlation is in agreement with Eq. (P^. 
In the thermal phase, where xq = in the large N limit, Eq. ( |100| ) reduces to 

= 1-12/- 1)91,, • 



in agreement with Eq. (^) and Eq. (185|). 



5.4 Statistical Barrier Penetration 

As we have seen above, the correlation length exhibits large peaks for the pump 
parameters 9q, Olf^_^_^ and/or The correlation grows as \/N at 6 = 9q. At 6lf^_^_^ and 
61i^ the large N dependence is, however, different. At these values of the pump parameter 
there is instead a competition between two neighbouring minima of the effective potential 
Vq{x) corresponding to, say, x = Xq and x = X2- A barrier, corresponding to a local 
maximum of Vo{x) at, say, x = Xi, then separates these two minima. Using the technique 
of [0 it can be shown that the peak in the correlation length close to = Olj,^^ is described 
by the expression 

27r 

7^ , , (105) 

[xi{l + nb) + bq{xi)] V-Vo"(xi) /(xo,Xi,X2) 

where we have defined 



/(xo,xi,X2) = v^K(^e"^[^o("i)-^o("°)l + V^K^^ • (106) 

We therefore conclude that 



7e^e^^^o ^ (107) 
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where AVq is the smallest potential barrier between the two competing minima of the 
effective potential Vo{x). This equation shows explicitly that the correlation length grows 
exponentially with N. Here we observe that Vo{6), and hence also AVq, is proportional 
to the combination a + nb{2a — 1) (see Eq. (|39D). Since the 7^-peaks grow exponentially 
with N they cannot be described by Eq. (^) with a finite power expansion of q{x). The 
first-order transitions at pump parameters corresponding to ^ = can be treated in a 
similar fashion. 



5.5 Discussion 



In Figs. |T0|, |TT] and [T^ we show the different analytical expressions for 7^ as given by 
Eqs. ([75D, ( ^41) and (|100|) . We have also plotted the exact correlation 7,^. For sufficiently 
large N, all the three expressions are good approximations to the correlation length around 
the first peak (see Figs. |10| and 0). We also observe that the mean field approximation 
is a better approximation than Eq. (0). Moreover, the mean field approximation as well 
as the the assumption 6^/x{9) + <^ 1 break down for pump parameters beyond 
the first peak of the correlation length. Hence, neither Eq. ( ^4|) nor Eq. (|100|) is a good 
approximation for 6''s above the first peak (see e.g. Fig. |ll|). Eq. ( |100| ), however, describes 
some qualitative features of the correlation length also for values of 6 above the first peak, 
even though Eq. ( p.OCI| ) then gives numerical values which are quantitatively wrong. We 
finally notice that all curves in Fig. O overlap in the large N limit. 



6 TRAPPING STATES 

The stationary photon probability distribution Eq. (||) has a special property when 
?7,fe = and Qm = 0. We then have that p„ = for all n > m. The cavity then cannot 
be pumped above m by photon emission from the pump atoms. The micromaser is then 
said to be in a trapping state f^, |16[-|1^. Actually, such trapping states have recently 



been observed in the stationary state of the micromaser system [^]. For a given A and 
k = 1,2, 3,..., we now define the function 

Vx + A^ 
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Figure 13: The solid line shows the order parameter (x) = {n)/N as a function 
of 9 when n?, = 0, a = 1, A = and = 100. As a guide to the eye we have 
also plotted the mean field solution (see Eq. ( p6|) and Eq. (p7|) ). The dashed-dotted 
curve is Xk{d) (see Eq. ( |109| )). For a given k > 1, this curve lies between the mean 
field solutions corresponding to the branches k and k — 1. 



or equivalently 



x,{e)=^-^-A^ . (109) 



For a given A^, trapping states then occur at the pump parameters 



OZk=(^k{Xtr) , (110) 

where Xtr = m/N and m = 0, 1, 2, 3, ... . 

The effect of trapping states on the order parameter (x) is illustrated in Fig. |T^ with 
an atomic flux parameter N = 100. We then clearly observe dips in the order parameter 
(x) due to the presence of trapping states. The observed structure of the dips in (x) can 
be explained as follows. For a given A; > 1, one can prove that the curve described by 
Eq. (|109|) lies between the mean fleld solution corresponding to branch k and k — 1 (c.f. 
Fig. |13D. Because of this fact, and since (n) < m, the trapping-dip at 6*^;. reaches down 
to (in fact, just below) the mean fleld curve corresponding to branch k — 1 (c.f. Fig. pIB]) . 
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In general, for a given k, trapping states such that the value of m/N is larger than 
the order parameter (x) have a minor effect on the order parameter. When 6 < 6^^ (or 
6fi) the presence of trapping states therefore have little effect on the order parameter (see 
e.g. Fig. |13D since the curve described by Eq. (|109D lies above the order parameter (x). 

It is easy to realize that trapping effects become less significant when the system 
is detuned. In this case we see from Eq. (^6]) and Eq. ( |109| ) that Xk{0) and the mean 
field solution is just reduced by the same amount A^. Numerical studies also show that 
the micromaser phase transitions occur essentially at values of the pump parameter 6 
almost independent of the value of the detuning < 1. Due to these facts the order 
parameter (x), Xk{0) in Eq. (|109|) as well as the mean field solution will intersect with 
a shifted horizontal ^-axe and hence trapping effects become less important. In passing 
we notice that, for a given value of the detuning A, trapping effects cannot occur for 



6 < = tt/vT+A^. This is due to the fact that the order parameter corresponding 
to the mean field solution (see Eq. and Eq. (0)) can never exceed unity. 

As the atomic flux increases, the dips in (x) due to trapping states become denser 
as a function of the pump parameter 6 [|T^. In the large N limit and for Uf, = 0, (x) 
therefore ceases to be an appropriate order parameter and the system appears to be 
frustrated. Since the observable (x) then varies rapidly it is natural to ask whether it is 
well defined at all. When the system is in a maser phase and when is sufficiently large, 
numerical studies of the standard deviation Ax = a/ (x^) — (x)^ as a function of 6 show 
that Ax is much smaller than the value of (x) itself for all possible values of 6, except for 
values of 9 near the phase transitions 9 = 9li^^^ and/or 9 = 91^^. Hence, the observable 
(x) is actually well defined except near these particular pump parameters. 

In Fig. |l^ we have studied the correlation length for different values of A^. The effect of 
trapping states becomes more pronounced when A^ increases. We observe that the peaks 
in the correlation length corresponding to trapping states get denser for large values of 
A^. Because of the same reasons as discussed earlier in this section, the effect of trapping 
states on the correlation length become less significant when the system is detuned. 

As mentioned in Section ^, the correlation length reaches a maximal value when the 
potential barrier between the two competing minima of the effective potential Vo(x) is at 
its lowest value. If the system is not in a trapping state, the probability distribution p„ 
then essentially consists of two Gaussian peaks with heights which are of the same order of 
magnitude. Furthermore, due to the pre-factor w{x) in p{x) (c.f. Eq. (|TB|)), a micromaser 
system which is close to a trapping state will also have a probability distribution which 
essentially consists of two Gaussian peaks of the same order of magnitude, even though the 
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Figure 14: The logarithm of 7^ as a function of 6 when = 25, = 50 and 
N = 100 and where = 0, a = 1 and A = 0. 

potential barrier between the two competing minima of the effective potential Vq{x) will 
then be smaller. Hence, the correlation length again reaches a large value. Exactly at a 
trapping state, however, the probability distribution consists essentially of one dominating 
Gaussian peak in the large limit. This means that we should expect the correlation 
length to become large close to a trapping state however not exactly at the occurence of 
such a state. Numerical studies are in accordance with these observations. As illustrated 
in Fig. |l5|, there is, furthermore, no visible difference between the ^-positions of the peaks 
of 7^ and the ^-positions of the dips of the (x) . This is not an obvious fact since ^ measures 
long-time features of the micromaser system but (x) is an instantaneous observable. As 
the number of branches increases, that is, the number of local maxima in p{x) increases, 
the connection between the heights of the peaks of ^ and the depths of the dips of (x) 
becomes increasingly more complex. 

As in Section |, let A„ denote the eigenvalues of the matrix L. Introducing the 
cumulative probability 

n-l 

Pn = J2P^ ' (111) 

it was then shown in Ref.0 that 
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Figure 15: The logarithm of 7^ and 15(x) as a function of 9 when the parameters 



are as in Fig. |13|. There is no visible difference between the ^-positions of the 
correlation trapping-peaks and the ^-positions of the order parameter trapping- 
dips. The dotted lines are the mean field solution Eq. ( p6|) and Eq. (p^) for branch 
k = 0, k = 1 and k = 2 for some limited interval of 6. 



( Pn{l-Pn) _ l-K/(l + n,r \ 

where ~ is an eigenvalue corresponding to the thermal equilibrium distribution. In 
the large N limit the left-hand side can be approximated by 7^ — 1. Hence, Eq. ( p.l2| ) 
reads 

V[(l + n,)n + iV6g„]p„ n ) ' ^ ' 

When the parameters are as in Fig. |T3| we find that, except for a small interval 2 < 
6' < 4, the sum-rule prediction Eq. ( |113| ) is in surprisingly good agreement with the exact 



correlation length, even at trapping states. We have no deep explanation of this numerical 
observation. 
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7 FINAL REMARKS 



We have studied the micromaser phase transitions at non-zero detuning and for pump 
atoms prepared in a diagonal statistical mixture. New novel features of the micromaser 
system then emerge as plateaux in the order parameter (x) as a function of the atomic 
density matrix as well as a twinkling behaviour at non-zero detuning. By introducing 
fluctuations in the pump-parameter 6 one can decrease the signal-to-noise ratio in observ- 
ables like V{+) and V{+, +,t) 0, Elsewhere we will investigate this feature of noise 

synchronization in more detail, also when the detuning A parameter is a random variable. 
It turns out that fluctuations in A also lead to features similar to fluctuations in 6, as 
was already anticipated in Ref. Fluctuations in the atomic density matrix parameter 
a lead, however, not to an output signal of the micromaser with less fluctuations. 
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